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I. INTRODUCTION 

It is commonly accepted that QCD is the correct the- 
ory for the strong interaction in the energy regime of 
interest for nuclear physics. At the same time, this is 
of limited practical value, because in this energy regime 
QCD needs to be solved non-perturbatively. Lattice sim- 
ulations, which arc in principle able to deal with this 
non-perturbative character are, for systems with A > 2 
nucleons, still in their beginning stages P, 0] ■ 

A possible way out of this dilemma is the applica- 
tion of chiral perturbation theory (ChPT) to nuclear sys- 
tems 0, HI' Making use of the spontaneously-broken 
chiral symmetry of QCD, one can formulate an effec- 
tive field theory (EFT) involving nucleons (W) and the 
Goldstone bosons related to chiral-symmetry breaking, 
the pions (it). ChPT is a powerful approach because 
it relates processes with different numbers of pions. In 
the purely-pionic and one-nucleon sectors, the Goldstone- 
boson character of pions guarantees that amplitudes can 
be expanded in powers of momenta 0. Q . In the few- 
nucleon sector, however, the existence of bound states 
clearly shows the non-perturbative nature of the prob- 
lem. 

Weinberg 0, @ recognized that this is caused by an 
infrared enhancement in the propagation of two or more 
nucleons. He suggested that the calculation of a generic 
nuclear amplitude should consist of two steps. In the 
first step, one defines the nuclear potential as the sum of 
"irreducible" sub-diagrams that do not contain purely- 
nucleonic intermediate states, and truncates the sum ac- 
cording to a simple extension of the standard ChPT 
power counting. In a second step, the potential is iterated 
to all orders, which can be done by using the Lippmann- 
Schwinger (LS) or Schrodinger equations. The potential 
includes pion exchanges and contact interactions, which 
represent the contributions of more massive degrees of 
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freedom. Assuming that contact interactions obey naive 
dimensional analysis, only a finite number of pion ex- 
changes and contact interactions contribute to the poten- 
tial at any given order: two non-derivative two-nucleon 
contact interactions and one-pion exchange (OPE) in 
leading order, and derivative contact interactions, two- 
pion exchange (TPE) and more-pion exchanges, and few- 
nucleon interactions in subleading orders. 

This power counting naturally explains that two- 
nucleon interactions are more important than three- 
nucleon interactions, etc. 9]. The resulting two-nucleon 
(AW) H3, E3 and three-nucleon (3W) [lj, El poten- 
tials provide a qua ntitative des crip tion of few-nucleon 
systems O EE EE E3, EE E3, EH H3, H3| In ad- 
dition, this approach matches well with the Nijmegen 
energy-dependent partial-wave analysis (PWA) of WW 
scattering data [24(. In this analysis, a separation of 
long- and short-distance physics is implemented by solv- 
ing the partial-wave Schrodinger equation with a long- 
range potential that consists of OPE and TPE (and the 
electromagnetic interaction), and a boundary condition 
with as many short-range parameters as are needed for 
an optimal description of the observables. The pion mass 
and OP E p arameters [25l |26| and even TPE parame- 
ters [23, |28| could be determined from the WW scattering 
data, in good agreement with values obtained from pion- 
nucleon scattering [23, yCJ |3l| . 

However, Weinberg's power counting has been criti- 
cized. As in any EFT, a regularization procedure is re- 
quired in order to separate high- and low-energy physics. 
Since this separation is arbitrary, a consistent power 
counting should provide sufficient counterterms at each 
order to absorb any cutoff dependence in the limit of large 
cutoffs. Because the solution of the LS equation is numer- 
ical in character, an explicit check of cutoff independence 
is challenging. This led Kaplan and coworkers [UIUOSI 
to examine a few of the diagrams contributing to the WW 
T matrix. They identified in two-loop diagrams ultravi- 
olet divergences proportional to the square of the pion 
mass and of the external momenta, which are present 
in leading order but cannot be absorbed by the avail- 
able counterterms. They concluded that pion exchange 
should not be fully iterated, but instead be treated in 
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finite order in perturbation theory. Quantitative calcu- 
lations at higher order showed, however, that this idea 
fails in some partial waves at momenta comparable to 
the pion mass |35j |. 

For smaller momenta, one can integrate out the pion 
and construct a "pionless" EFT, which is very successful 
within its limited range Q|. Nevertheless, a lot of interest- 
ing nuclear physics is thought to take place at momenta 
of the order of the pion mass. (The Fermi momentum of 
isospin-symmetric nuclear matter, for example, is about 
300 MeV.) It seems unavoidable that in this larger mo- 
mentum range pion exchange has to be iterated. 

It is now well known that the renormalization of an 
EFT is not necessarily the same as that of its perturba- 
tive series. This is seen clearly in the three-body problem 
in the pionless EFT j3|||32], OHI . The origin of this feature 
lies in the renormalization of singular potentials [39l Eof . 
In the specific case of OPE, the singularity is the 1/r 3 
behavior of the tensor force in spin-triplet channels. It 
has been found that the cutoff dependence of an uncou- 
pled 1 / r n interaction in the S wave can be absorbed into 
one counterterm |4(ij |. 

The renormalization of OPE in lowest waves was re- 
examined from the non-perturbative viewpoint in Refs. 
[4ll I42] ]. The problem with the ultraviolet divergence 
proportional to the pion mass squared in the 1 So chan- 
nel persists in this context O n the other hand, the 
divergence associated with momenta, present in the 3 Si- 
3 Di coupled channel, can be absorbed into the available 
leading-order counterterm [ill . With a further ex- 
pansion around the chiral limit, Weinberg's power count- 
ing seems to be consistent in a non-perturbative calcula- 
tion of these waves 0> D3 13 • (F° r a different conclu- 
sion, see Ref. |45|.) 

OPE contributes, however, also in higher partial waves. 
The naive power counting does not predict leading-order 
counterterms in these partial waves. However, the sin- 
gularity of the tensor interaction exists in all the spin- 
triplet channels. In fact, it has been argued that for an 
uncoupled singular interaction boundary conditions need 
to be fixed in all waves where the potential is attrac- 
tive 0. Therefore, cutoff dependence can be expected 
in some spin-triplet channels if there are no correspond- 
ing counterterms, posing a significant difficulty for Wein- 
berg's power counting. 

Another important renormalization issue concerns few- 
nucleon forces. In the pionless theory it has been shown 
that consistent renormalization requires a 3N force in 
leading order [3^]. This result is not necessarily in con- 
tradiction with Weinberg's power counting in the "pion- 
ful" EFT, because the 3 AT force in the pionless theory in- 
cludes contributions that are iterations of the AW force 
with intermediate-state nucleons of momentum 0(m n ) 
in the pionful theory. The two EFTs have AW interac- 
tions with different ultraviolet behaviors. Whether OPE 
sufficiently softens the asymptotic behavior of the 3N LS 
equation is an issue that remains unresolved. 

In practice, the renormalization issue has been 



sidestepped by choosing rather low cutoffs to regular- 
ize the LS equation and by varying the cutoffs only in a 
very limited range [13, [IE EE 113, 113, . 

Cutoff dependence has generally been observed in higher 
partial waves; see, for example, the discussion in Ref. 
[43 regarding higher orders in the Weinberg expansion. 
For relatively small cutoff variations, it has been noticed 
that the resulting variations in the phase shifts decreases 
with increasing order. It has then been assumed that the 
observed cutoff dependence is of the order of the error 
expected from the truncation of the expansion. In fact, 
it has been argued that the EFT involving nucleons and 
pions necessarily involves a mild cutoff dependence and 
that cutoff values exist that are optimal for the conver- 
gence of the expansion 0> E3 • 

An in-depth study of the cutoff dependence in higher 
AW partial waves and in the 3A~ system still needs to be 
performed. This study is the aim of this work. We seek 
to quantify the cutoff dependence in lowest order and, 
if possible, to identify ranges of cutoffs in which only 
small variations of observables occur. We then discuss 
how to absorb the cutoff dependence in a finite num- 
ber of counterterms. We first consider AW scattering, in 
which case we compare our results with the phase shifts 
and mixing angles from the energy-dependent Nijmegen 
PWA93 [24( , which provide an optimal representation of 
the AW database (The fact that this PWA93 did not 
yet include TPE [27I |28| does not affect our investiga- 
tion of OPE.) We then extend our analysis to the 3AT 
bound state. We will restrict ourselves to total AW an- 
gular momentum j < 4, which is sufficient to study the 
3N binding energy. As we are going to show, our results 
have significant implications: we present here a modifica- 
tion of power counting that is consistent with all known 
results, and could become a new basis to organize nuclear 
interactions in the EFT with pions. 

Section [H] describes the interaction and our approach 
to regularize and solve the AW Lippmann-Schwinger 
equation. In Section II I II we identify problematic par- 
tial waves, explicitly show their cutoff dependence, and 
present counterterms that generate cutoff-independent 
phase shifts in reasonable agreement with the PWA. Sec- 
tion ^3 is devoted to the 3N system. The implications 
of our findings to power counting in nuclear ChPT are 
analyzed in Section Finally, our conclusions and an 
outlook are given in Section IVT1 



II. REGULARIZATION OF THE 
LIPPMANN-SCHWINGER EQUATION 

We first consider AW scattering in the center-of-mass 
frame. We denote by /i the reduced mass (mjv = 2fi is 
the nucleon mass), by E the energy, and by p and p' 
the relative momenta before and after interaction; the 
momentum transfer is q — p — p' . The relative distance 
between the two nucleons is r. The standard Pauli ma- 
trices in spin and isospin space are denoted by <r, and t, , 
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respectively. 

With a standard normalization for plane waves, the LS 
equation for the T matrix reads, in momentum space, 



T(p',p,E)=V(p',p) 
' d 3 p"V(p\p")- 
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with V the potential. The OPE potential is 
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where is the pion mass. In lowest order the strength 
of OPE is completely determined by the axial-coupling 
constant qa — L26 and the pion-decay constant f n = 
92 A MeV. 

In addition to pion exchanges, the EFT contains short- 
range interactions that represent high-energy degrees of 
freedom that have been integrated out. The simplest are 
two contact interactions 
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(c p s + a P t ) . 
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where we used the projectors onto spin-triplet and spin- 
singlet states, Pt and P s . The two strength parameters c s 
and Ct need to be determined from NN scattering data, 
for instance from the scattering lengths in the 1 So and 
3 Si channels. It is possible to write 
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where the parameters Cq and D 2 are independent of the 
quark masses. 

For the numerical solution of the LS equation, we 
need to introduce a regulator f(p',p) that effectively 
cuts momenta at a cutoff A. The regularization proce- 
dure is an arbitrary splitting of short-range physics into 
the high-momentum region of loops and contact inter- 
actions. Low-energy physics should, of course, be inde- 
pendent of the choice of regulator (renormalization-group 
invariance), once the dependence of contact parameters 
on the cutoff is taken into account. It is convenient for 
the partial-wave decomposition to perform the regular- 
ization using momentum cutoff functions depending on p 
and p' rather than on q. Here we use 



f(p',p) = e 



-(p 4 +p' 4 )/A 4 



(5) 



This leads to nonlocal interactions in configuration space. 
However, because the regulator only depends on the mag- 
nitude of the relative momenta, it does not influence the 
partial-wave decomposition. This guarantees that con- 
tact interactions act in specific partial waves, indepen- 
dent of A. In particular, it implies that V c only acts in 
the two S waves. 

For the following discussion, it is useful to look also at 
the configuration space expression for OPE, 



Ti ■ r 2 [T{r) S12 + Y(r) d\ ■ a 2 ] , 

(6) 
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TABLE I: Matrix elements of the operator n -Ta 5*i2 for spin- 
triplet channels with total angular momentum j. The matrix 
elements depend on the isospin t, and on the incoming and 
outgoing angular momentum / and 
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r(r) = 
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and the tensor operator is 

#12 = 3(cti • f)(a 2 
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The tensor force T(r) of OPE contains a singular inter- 
action ~ 1/r 3 that acts in the spin-triplet waves; the 
tensor force is zero in the spin-singlet channels. Using 
the partial-wave matrix elements given in Table^one can 
identify whether the tensor force is attractive or repulsive 
in specific partial waves. This we will require below. 



III. NUCLEON-NUCLEON PHASE SHIFTS 

Our aim is to study the dependence of observables on 
the chosen value for the cutoff A. We have performed a 
partial-wave decomposition of the interaction described 
in the previous section and then solved the LS equation 
and extracted phase shifts. The explicit expressions are 
summarized in Appendix lAl We study the cutoff depen- 
dence of the phase shifts in leading order (LO), or O(Q ), 
in Weinberg's power counting. We consider A in a wide 
range, between 2 fm _1 and 20 fm _1 . 

We start with the S-wave channels, which were pre- 
viously examined in Refs. 0, E3, E3> Lul with different 
regularizations. We fit c s and Ct to the 1 Sq and 3 Si phase 
shifts at 10 MeV and we confirm the cutoff independence 
found in Refs. |4lll42|. as can been seen in Figs. ^ and [21 
In Fig.^we show the running of c s with the cutoff A, and 
the resulting cutoff dependence of the 1 So phase shifts 
at various laboratory energies. In Fig. [21 we show the 
corresponding results for ct, and the 3 Si and 3 Di phase 
shifts and the mixing angle E\. One sees that the cutoff 
dependence of the phase shifts is small for A>5 fm _1 , 
but it increases for increasing energy, as expected in an 
EFT. It is interesting to note that Ct(A) displays a nice 
limit-cycle-like behavior, similar to the 3A^ force in the 
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3iV problem in pionless EFT [33] , which is solved using a 
regulator similar to ours. Since the counterterm strength 
behaves differently in Ref. 0] , where a coordinate-space 
regulator was employed instead, we conclude that this 
behavior is regulator dependent. This is in line with a 
similar recent finding for a central potential |HcLl5lj . Note 
that the running found here is similar to that observed 
in Ref. |52| for a counterterm in a different channel at 
higher orders but with the same regularization of the LS 
equation. This suggests that the form of the running is 
perhaps more influenced by the regulator than the spe- 
cific form of the singularity of the interaction. 

The resulting phase shifts as function of the laboratory 
energy are shown in Fig. |3J In the channel, we re- 
cover the known strong deviation of the LO result from 
the PWA. This is related to the relati vely -large effective- 
range parameter in this partial wave |42l l53| . which can- 
not be reproduced without a two-derivative contact inter- 
action. This problem is solved once the latter is included 
in subleading order (see, e.g., Ref. ^3)- ^ n the coupled 
3 Si- 3 Di channels, we find an encouraging agreement with 
the PWA. The mixing angle e\ is underpredicted, when 
one goes to the limit of high A. In this limit the agree- 
ment with the PWA has, however, improved c omp ared 
to the choice A ~ 3 fm -1 used in the literature |16| . 

Despite these positive results, there are potential prob- 
lems in other waves. The explicit expression in Eq. 10 for 
Vitt in configuration space suggests that, because of the 
singularity of T(r), we can expect a number of bound 
states — infinite in the limit A — > oo — in all channels 
where the tensor force is attractive. A consequence would 
be cutoff dependence in these waves. We will, therefore, 
study higher partial waves, starting from P waves, since 
these, according to Weinberg's power counting, should 
not require counterterms in lowest order. 

We start with spin-singlet channels, where the tensor 
interaction is zero and the potential is non-singular. The 
results for the 1 Pi, X D 2 , 1 F3, and : G4 phase shifts as 
function of A at four different lab energies are shown in 
Fig. 0] It is seen that indeed the dependence on A be- 
comes smaller with increasing A. Even for a rather high 
energy of 190 MeV, we find in all cases only negligible 
changes in the phase shifts for A>5 fin -1 . This sup- 
ports the claim that in these channels no inconsistency 
in Weinberg's power counting exists. Fig. [5] compares 
the resulting phase shifts to the PWA. For the P and D 
waves the agreement is good below 30 MeV. Above that 
energy significant higher-order contributions are neces- 
sary to improve agreement with the PWA. For the F and 
G waves, where contact interactions are expected at even 
higher orders, the agreement is much better for energies 
up to 100 MeV. We note that for the energies below 30- 
50 MeV, where all singlet phase shifts are described well 
by the lowest-order predictions, the cutoff dependence is 
already negligible for A > 2 frn -1 . 

In the next step, we look at the triplet channels where 
the tensor force is repulsive. The 3 Pi and 3 F3 partial 
waves belong to this class. The A dependence of these 



phase shifts is shown in Fig. and, again, we obtain A 
independent results for A>5 frn -1 even for energies as 
high as 190 MeV. The comparison with the data for these 
cases is shown in Fig. [7] The 3 Pi result is in much better 
agreement with the PWA than the corresponding result 
of the 1 Pi singlet channel. The F wave has a compa- 
rable accuracy in the triplet and singlet cases. These 
results confirm that Weinberg's power counting is, again, 
consistent in channels without an attractive singular in- 
teraction. The predictions agree well with the low-energy 
data in all these cases. 

Next, we look at the triplet channels in which an at- 
tractive singular tensor force acts, and therefore the un- 
regulated problem is not well defined. At a finite cut- 
off bound states occur in various waves, and increasing 
A generates more and more of these bound states. For 
A between 2 and 20 frn -1 , we find bound states in the 
3 P and 3 D 2 channels, see Fig. |S] In the 3 P 2 - 3 F 2 par- 
tial waves, we find a bound state for cutoffs just above 
A = 20 frn -1 . In higher partial waves the interaction 
appears to be screened enough by the centrifugal barrier 
that no bound state occurs in our cutoff range, although 
bound states should appear at sufficiently-higher cutoffs. 

Even though the binding energies increase rapidly with 
the cutoff, the bound states appear at zero energy. It 
is to be expected that in the cutoff regions where new 
bound states appear the variations of the phase shifts 
are strong. This is explicitly shown in Fig. [3] for the 
phase shifts at 10 MeV and 50 MeV in various attractive 
channels. Clearly, an inconsistency in Weinberg's power 
counting shows up in these channels, because there are no 
counterterms available to remove the cutoff dependence 
of the observables. 

It is interesting to note that in some cases we can iden- 
tify cutoff regions in which the results are stable and all 
the bound states are deep. Some clear plateau regions oc- 
cur, especially at the lower energies. At 50 MeV, the 3 D 2 
phase shift in the plateau region is ~ 9°, which agrees 
very well with the Nijmegen PWA (8.97°). The corre- 
sponding 3 D3 phase shift, however, is too small. The 
situation is even worse in the P waves. At the same en- 
ergy, the cutoff dependence of the 3 P and 3 P 2 phase 
shifts remains visible in any region of As. 

This cutoff dependence is related to the singularity of 
the interaction. It is known that an attractive singular 
central potential requires a boundary condition in each 
partial wave |4^. Therefore, we propose to add to each 
of the problematic triplet channels a counterterm and 
fit it to experiment, say to the phase shift at a certain 
(low) energy. We then show that the cutoff dependence 
indeed vanishes also for other energies. In the following, 
we will illustrate this explicitly for the 3 Po, 3 D 2 , and 3 P 2 - 
3 F 2 partial waves, which are the most problematic cases, 
because bound states exist or are close to appearing in 
the cutoff range that we examined. Our results extend 
the work of Refs. [U |H |H Q to channels beyond S 
waves (and to our choice of regulator). 

As argued, we add contact interactions in the 3 Pq (i = 
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FIG. 1: Fit result for the counterterm c s as a function of the cutoff, and the resulting cutoff dependence of the So phase shifts 
at laboratory energies of 10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 190 MeV (dash-dotted line). 



120 
100 
W> 80 

D 

5? 60 
40 
20 



I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 



'I i I i I i I i I i I i I 



0.0 
-5.0 



■— , -10.0 
-S -15.0 



to 



■20.0 
■25.0 
■30.0 





2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 

A [fm" 1 ] A [fm" 1 ] 

FIG. 2: Fit result for the counterterm c t as a function of the cutoff, and the resulting cutoff dependence of the 3 Si- 3 Di phase 
shifts and the mixing angle ei at laboratory energies of 10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 
190 MeV (dash-dotted line). 



1) and 3 P2- 3 F 2 channels (i = 2) of the form 
1 c 
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which in Weinberg's power counting appear only at next- 
to-leading (NLO) order, or 0(Q 2 ). The first D-wave 
counterterms are supposed to be of even higher order: 
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Fig. 1 11)1 shows our result for the 3 Po partial wave. The 
value of C\ was determined by a fit of the phase shift 
for a laboratory energy of 50 MeV. Since the size of the 
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FIG. 3: Comparison of the 1 So and 3 Si- 3 Di phase shifts and the mixing angle ei (as function of the laboratory energy) in 
lowest order for A = 20 fm _1 (solid line) to the Nijmegen PWA (dashed line). 



counterterm is not bounded, we varied this constant by 
orders of magnitude, but could not find any further so- 
lution that describes the phase shifts equally well. The 
cutoff dependence of c\ exhibits a nice limit-cycle-like be- 
havior, similar to that of ct- Fig. EH also demonstrates 
that the resulting phase shifts at other energies are cutoff 
independent for A>8 fin -1 . Figs. II II and II 21 summarize 
the analogous results for the 3 P2- 3 F2 and 3 D2 partial 
waves, respectively. The fits were performed using the 
3 P2 phase shift at 50 MeV and the 3 D 2 phase shift at 
100 MeV. We confirm the cutoff independence (for large 
A) in all phase shifts and mixing parameters. 

An alternative to absorbing the cutoff dependence in 
the various P waves individually would be to employ one 
counterterm with tensor structure. Unfortunately, we 
have not been able to implement this idea without intro- 
ducing cutoff dependence in the 3 Pi wave. 

After removing the cutoff dependence by adding appro- 
priate counterterms, we still find spurious bound states 
in the 3 Po, 3 D2, and also the 3 Si- 3 Di channels. However, 
the cutoff dependence of the binding energies is now com- 
pletely different, as shown in Fig. ^5] As desired, only 
3 Si- 3 Di has a shallow bound state, the deuteron, which 
is cutoff independent over almost the entire A range; the 
deuteron binding energy is predicted to be 1.92 MeV in 
this LO calculation. The bound states in the other chan- 
nels are all very deep. A new bound state appears with 



infinite binding energy around the cutoff at which the cor- 
responding counterterm is singular, and then approaches 
a constant, large binding energy for increasing A. These 
bound states are beyond the range of the EFT, and they 
are irrelevant for the low-energy physics. 

With the added counterterms, we obtain a very decent 
description of the phase shifts. Fig. 1141 shows that our 
3 Po result follows the energy-dependence of the Nijme- 
gen PWA remarkably well. Obviously, the addition of 
the counterterm is here supported by the experimental 
data, fn the coupled 3 P2- 3 F2 channels the agreement 
with the PWA below 50 MeV is still satisfactory. We 
emphasize that the 3 F2 phase and the mixing parameter 
£2 are predictions. Choosing a high cutoff A clearly does 
not compromise the description of these observables. 

For the 3 D2 phase, see Fig. we find again a good 
agreement with the PWA. Here, we also included the 
prediction based on a calculation without counterterm, 
for A = 8.0 fm _1 in the plateau region of Fig. [5] Fol- 
low energies below 50 MeV, the results are comparable. 
The deviations from the PWA become significant toward 
higher energies, where the plateau seen in Figgis more 
and more tilted. For these higher energies, the counter- 
term again improves the predictions. 

Our overview is completed in Figs. IT51 and lItjl with the 
3 D 3 - 3 G 3 , 3 F 4 - 3 H 4 , and 3 G4 channels. In these partial 
waves there is a relatively small cutoff dependence in the 
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FIG. 4: Cutoff dependence of the singlet phase shifts for various partial waves. Results are given for laboratory energies of 
10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 190 MeV (dash-dotted line). 



A range that we studied (although presumably cutoff de- 
pendence will become significant at cutoffs high enough 
to bring in spurious bound states) . In all cases the agree- 
ment with the PWA is improved when we increase the 
cutoff from the traditional values around 2.5 fm -1 [l^ 
to our higher values. This is especially true for the 3 D3 
partial wave, which, for our higher cutoffs, becomes at- 
tractive for higher energies. 

After these encouraging results, we examine the 3N 
bound state in the next section. 



IV. THREE-NUCLEON BOUND STATE 

The power of EFT comes to bear when more nucle- 
ons are considered. The 3A" system is the first extension 
to few-nucleon systems to consider. According to Wein- 
berg's power counting, 3 AT forces should be subleading in 
the pionful theory, even though there is a leading-order 
3N force in the pionless EFT [2a|. Calculations of the 
triton properties using traditionally-low values for A have 
been published in Refs. Q for NLO and N 2 LO 

in Weinberg's power counting, with and without 3N in- 
teractions. So far, a LO calculation was omitted, because 
of the unsatisfactory description of the AW phase shifts, 
especially So- Here, our main goal is to assess the or- 



dering of 3A~ counterterms. For this purpose we do not 
require a high-quality description of the NN phase shifts, 
but instead we need to study the dependence of the triton 
binding energy on a larger range of cutoffs. 

The 1 So problem cannot be addressed in this work, 
but a calculation of the triton binding energy (E t ) is of 
interest in order to compare the renormalization of the 
3N system in the pionful theory with that in the pion- 
less EFT. We will demonstrate that no 3N counterterms 
are necessary to ensure cutoff independence in LO in the 
pionful theory, once the LO calculation has been modi- 
fied according to the previous section. This result will be 
important in future studies of the convergence of the chi- 
ral expansion in few-nucleon systems. The appearance 
of 3N interactions at N 2 LO, the parameters of which are 
generally adjusted to the experimental value of E t , makes 
predictions in LO and NLO especially important in this 
respect. We defer a study of subleading orders and Nd 
scattering states to a later publication. 

We calculate E t by solving the Faddeev equations for 
the 3N system. The techniques were recently described, 
e.g., in Ref. [54|. Here we have to deal with the ad- 
ditional complication that the NN interaction supports 
deep, spurious bound states in various partial waves, 
which we remove as described in Appendix [5] We con- 
firmed the accuracy of this prescription by comparing 
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FIG. 5: Comparison of the LO singlet phase shifts (as function of the laboratory energy) for A = 20 fm 1 (solid line) to the 
Nijmegen PWA (dashed line). 
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FIG. 6: Cutoff dependence of the repulsive triplet phase shifts for two partial waves. Results are given for laboratory energies 
of 10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 190 MeV (dash-dotted line). 



our results for the energy to the expectation value of the 
Hamiltonian using the unaltered interaction. Both values 
agreed within several keV. 

The cutoff dependence of E t is shown in Fig. ~ 7| We 
see a plateau region starting around A = 8 fm -1 . To 
extract the converged result from this calculation and to 
confirm that only terms of order 0(Q 2 /A 2 ) and higher 
are missing, we fitted the function E(x) = Eq (1 + 



(C/A) x ) to our numerical results. We obtained the con- 
verged binding energy E Q = —3.6 MeV, C — 2.54 fm -1 , 
and x ~ 1.8. The exponent is in reasonable agreement 
with the expectation x — 2. The quality of this approxi- 
mation to the cutoff dependence can be observed in the 
Lepage plot j - J of Fig. ~ ~ where we show E — Eq versus 
A on a double-logarithmic scale. The data follows the 
fitted results nicely. For high values of A the slope seems 
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FIG. 7: Comparison of LO repulsive triplet phase shifts (as function of the laboratory energy) for A = 20 fm 1 (solid line) to 
the Nijmegen PWA (dashed line). 
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FIG. 8: Binding energies of the spurious bound states in selected attractive triplet channels. 



to change slightly. This can be a numerical artifact, be- 
cause E—Eq is rather small in this range, which increases 
the relative, numerical uncertainty. 

A cutoff dependence as it occurs in the 3N system 
within the pionless EFT [3|| does not appear here. We 
stress that, in contrast, for the same range of A we found 
considerable variations in several NN phase shifts before 
adding counterterms. This gives us confidence that our 
cutoff range is large enough to draw conclusions about 
the order of the counterterms. In particular, we con- 
clude that the finite range of OPE prevents the Thomas 
collapse of the 3 N bound state |5^| . 

Our LO results imply a rather-large underbinding of 
the triton in the limit of high A. In the region around 
A = 2.5 fin , our prediction is in much better agree- 
ment with the experimental value E t = —8.48 MeV. For 
those cutoffs the Weinberg NLO predictions are very sim- 
ilar |2lJ . This suggests that the NLO result might be less 
cutoff dependent than the LO result, and that the plateau 
region starts at lower values of A for higher orders. Our 
result is probably influenced significantly by the unre- 
alistic description of the AW phase shift, which will 
improve in a NLO calculation. Binding energies are more 



sensitive to higher orders, because in a theory with pions 
the potential energy is to a large extent canceled by the 
kinetic energy. In view of this, we are not very concerned 
about our LO plateau value. It will be interesting to ana- 
lyze the NLO results in a similar way. If our expectations 
are verified, they would reconcile our observations with 
the good results obtained with a fixed, low cutoff in Refs. 

Mm 

V. LESSONS FOR POWER COUNTING 

It has been realized for some time now that EFT power 
counting is more complicated for non-perturbative than 
for perturbative processes. In particular, one has to con- 
sider the effects of infrared enhancements in the running 
of counterterms, which can invalidate naive dimensional 
analysis. In this section, we discuss some of the impli- 
cations of our findings to power counting in the pionful 
EFT. 

The existence of shallow (real and virtual) bound 
states in both NN S waves is a clear sign that non- 
perturbative physics is taking place, in contrast to the 
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FIG. 9: Cutoff dependence of phase shifts in attractive triplet channels at laboratory energies of 10 MeV (solid line), 50 MeV 
(dashed line) , and 100 MeV (dotted line) . 
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FIG. 10: Fit result for the counterterm ci as a function of the cutoff, and the resulting cutoff dependence of the 3 Po phase 
shift at laboratory energies of 10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 190 MeV (dash-dotted 
line) . 



situation in ordinary ChPT. We can describe this in the 
same language used to discuss power counting in ChPT 
we represent typical nucleon momenta by Q and 
the characteristic scale of QCD in the hadronic phase by 
Mqcd- The effect of iterating an interaction in the ker- 
nel of the T matrix is twofold. First, one has an extra 
three-dimensional momentum integral and an extra NN 



Schrodinger propagator. Second, one has an extra factor 
of the potential. After the cutoff dependence is removed 
by renormalization, the contribution to the NN T ma- 
trix from an NN intermediate state is expected to be 
C(mAr(5/47r). This is a factor tum/Q ^> 1 larger than 
in analogous states in ordinary ChPT, and it is due to 
the small energy of intermediate states containing nuclc- 
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FIG. 11: Fit result for the counterterm C2 as a function of the cutoff, and the resulting cutoff dependence of the 3 P2- 3 F2 phase 
shifts and the mixing angle £2 at laboratory energies of 10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 
190 MeV (dash-dotted line). 
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FIG. 12: Fit result for the counterterm Cd as a function of the cutoff, and the resulting cutoff dependence of the 3 T>2 phase 
shift at laboratory energies of 10 MeV (solid line), 50 MeV (dashed line), 100 MeV (dotted line), and 190 MeV (dash-dotted 
line) . 



ons only. If the interaction has an effective strength via 
connecting waves of orbital angular momenta I and I', 
one iteration then roughly brings a dimensionless factor 
O^nQvi'i/^tt). For Q > iir/niNVin this factor is > 1 
and the interaction has to be iterated to all orders, po- 
tentially leading to bound states. 



The pionful theory is relevant for momenta Q > m n . In 
this case, we estimate the effect of OPE as vin <~ ctin/f%, 
where am is a dimensionless angular-momentum factor. 
Therefore, we can expect OPE to be non-pcrturbative at 
Q > Anf^/niNai'i- Ignoring the angular-momentum fac- 
tor, we get Q > 100 MeV for the momentum where pions 
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FIG. 13: Binding energies of bound states found in various 
partial waves as function of the cutoff. The shallow bound 
state in the 3 Si- 3 Di coupled channels corresponds to the 
deuteron, while the other, deep, bound states are outside the 
range of applicability of the EFT. 



need to be iterated. This is in agreement with what is ob- 
served in a perturbative calculation of low waves, where 
twice-iterated pion exchange seems to overcome one-pion 
exchange in various channels for momenta around 100 
MeV I23. 

In spin-singlet channels, OPE goes as 1 — m^/q 2 + . . . 
at high momentum. When iterated, the first term by 
itself introduces cutoff dependence in the S wave only, 
which can be removed by a chirally-symmetric counter- 
term Co- The interference between the iteration of this 
counterterm and the second term in OPE generates fur- 
ther cutoff dependence in the wave, which in turn 
can be removed [32I |42| by a chiral-breaking counter- 
term m 2 D2, cf. Eq. 0. This counterterm is enhanced 
with respect to naive dimensional analysis, and should 
be promoted to LO if pion-mass effects are kept in LO, 
as it seems most efficient. In this paper we have checked 
explicitly that OPE is well behaved in the higher spin- 



singlet waves. (Higher spin-singlet waves have also been 
recently discussed in Ref. |5q.1 

In spin-triplet channels, the situation is complicated 
by the tensor operator, which retains angular depen- 
dence even asymptotically. As shown in Refs. 0, 
l42| and confirmed here, the cutoff dependence intro- 
duced by iteration in the 3 Si- 3 Di coupled channel can 
be dealt with a single chirally-symmetric, momentum- 
independent counterterm, as prescribed by Weinberg's 
power counting. The deuteron and 1 So virtual state 
have, for the observed value of the quark masses, bind- 
ing momenta somewhat smaller than our crude estimate 
Q ~ 4irf 2 /mpf, indicating an amount of fine tuning. 
However, if one varies the pion mass the momentum 
scales for the bound states acquire more natural values 

mm mm. 

Yet, as we have shown above, iterated OPE produces 
spurious bound states and cutoff dependence in all waves 
where the tensor force is attractive. This undesired fea- 
ture can only be removed by additional counterterms in 
the corresponding waves. Since Weinberg's power count- 
ing only prescribes counterterms in the S waves, our re- 
sults clearly upset Weinberg's power counting. 

Weinberg's power counting was based on naive di- 
mensional analysis. Contact interactions are necessary 
to remove divergences from loops that do not involve 
purely- nuclconic intermediate states. These loops are 
not infrared enhanced, and it is reasonable to assume 
that they scale — as in ordinary ChPT — with powers of 
(Q/47r/ w ) 2 . This has in fact been confirmed by explicit 
calculation E3, HJ IP- Implicit in Weinberg's power 
counting is the assumption that loops in the LS equation 
do not bring significant new cutoff dependence. The pa- 
rameters of contact interactions with derivatives or pow- 
ers of the pion mass would thus be suppressed by powers 
of a large mass scale, Mqcd, and the effects of derivatives 
would scale as Q/Mqcd- 

However, we now see that Weinberg's implicit assump- 
tion is not correct. The short-range parameters needed 
to renormalize iterated OPE do not obey naive dimen- 
sional analysis because their renormalization-group run- 
ning is enhanced in the infrared. These counterterms are 
driven by pion parameters, and the effects of derivatives 
can scale as Q/f^. (If we take rnjv ~ 47r/ 7r , there is no 
other dimensionful parameter at LO than f w .) Taken at 
face value, this implies that all these counterterms must 
be considered in leading order. This is not a complete 
disaster, as there is still some predictive power left, for 
example in the energy dependence of each attractive par- 
tial wave, and in the repulsive waves. However, it would 
put few-nucleon observables that include significant con- 
tributions from many partial waves, such as the triton 
binding energy, out of reach. 

In the following, we want to motivate that this com- 
plication can be avoided in higher partial waves, because 
these are still perturbative. We therefore consider ain, 
which involves a kinematic suppression that accounts for 
the repulsive effect of the centrifugal barrier. The ap- 



13 




propriate counterterms will make OPE well defined, by 
selecting in the region r ~ 1//tt the correct combination 
of solutions for the long-range potential 0, 03 • There- 
fore, the kinematic suppression can be estimated as for 
a regular potential. In the case of a central potential, 
it can be shown that for I 3> Qd, where d is the 
range of the interaction, the l-wave phase shift is given 
by tan St ~ (Qd/(l + l/2)) 2l+1 < 1 (barring the excep- 
tional case of a fine tunning that generates a T-matrix 
pole near the origin of the complex momentum plane). 
This is consistent with the expectation that in the clas- 
sical limit there is little scattering when the impact pa- 
rameter l/Q is much larger than the range d. The ratio 
of the T matrix, and thus the potential, between / + 1 
and I is 0(Q/lm 7T ) 2 , for large I. For Q ~ m„, we are led 
to a>tj — 0(1/ ll 2 ). In the case of the tensor force, there 
must be two different factors 1//! and 1//'!, in addition 
to the elements in Table H] which do depend on I and I', 
but approach constant values for large 1,1'. Although at 
small I and /' ayi can be very complicated, we expect 
that ain = 0(l/l'W.) for large I, I'. For this argumenta- 
tion we assumed that the pion mass is finite. We note 
that it does not apply in the chiral limit. 

This qualitative argument suggests that the effects 
of the corresponding higher-derivative counterterms are 
suppressed by a large (for large I) scale lf n . Obviously, 



there might be other dimensionless factors that we miss 
here, but the fact that factors of I suppress OPE and its 
required counterterms in high-Z waves must hold. There 
are several implications of this new counting that seem 
to be supported by existing results. 

The kinematic suppression of higher waves makes the 
strength of OPE weaker with increasing In high waves, 
OPE and probably all subleading interactions can be 
treated in perturbation theory even for momenta of the 
order of . We therefore find explicitly that high partial 
waves can be treated perturbatively. The perturbative 
nature of OPE in high waves is part of nuclear folklore, 
and has been checked in EFT explicitly E01 EH ■ 

For sufficiently-large I, the suppression factor in 
counterterms becomes dominated by Mqcd (rather than 
lf n ), representing omitted QCD degrees of freedom, and 
the size of the counterterms is that assumed in Wein- 
berg's power counting. 

On the other hand, for a finite number of low partial 
waves, e.g. 3 Pq and 3 P2- 3 F2, we find that perturbation 
theory is not sufficient for momenta of the order of m^. 
This is caused by the lack of enough suppression from I, 
and by unnaturally-large \aw \ in these cases. Resumma- 
tion is necessary and needs to be performed numerically. 
Our numerical analysis showed that the cutoff depen- 
dence can be absorbed by one counterterm per partial 
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FIG. 15: Comparison of attractive triplet phase shifts (as function of the laboratory energy) for A = 20 fm 1 (solid line) to 
the Nijmegen PWA (dashed line). For 3 T>2 also the result for A = 8.0 fm -1 without counterterm is given (dotted line). 



wave. The favorable agreement of our LO calculation 
with the data indicates that we did not introduce addi- 
tional inconsistencies and that the resummation includes 
the most important diagrams. 

We conjecture that this mixture of perturbative treat- 
ment of higher partial waves, resummation of lower par- 
tial waves, and promotion of a finite number of counter- 
terms is the most consistent approach to ChPT for nu- 
clear systems. In subleading orders, it would natu- 
rally suggest a perturbative treatment of the sublead- 
ing interactions, if no unnaturally-large aws are present. 
The NLO interactions consist in principle of TPE and 
counterterms with two more derivatives than LO. Sub- 
sequent orders are constructed by the inclusion of suc- 
cessive powers of Q/Mqcd- The most-effective organi- 
zational scheme for subleading interactions probably re- 
lies on taking into account an explicit delta-isobar field 
[Iall4llMl6lLl6^ |. The correctness of our modified power 
counting needs, of course, to be checked in future studies 
of higher orders. 

On the other hand, for practical reasons, it might 
be convenient to perform the resummation in all par- 
tial waves that are taken into account. Iterating some- 
thing small causes only a small error, so one might de- 
cide to iterate OPE in all waves, as done automatically 
when solving the LS for OPE. This is again part of nu- 
clear folklore. If we do this without the corresponding 



counterterms, however, cutoff dependence is introduced. 
By increasing the cutoff, the iteration of OPE can be 
made arbitrarily large, and at some point bound states 
appear. The kinematic suppression suggests that the cut- 
offs for which bound states appear increase with I, which 
is consistent with what we observe in the cutoff window 
we studied. Existing calculations based on Weinberg's 
counting should exhibit some of the problems we 
point out here. In particular, as the cutoff is increased, 
partial waves without the required counterterms become 
unwieldy. 

In that case, one should consider cutoffs in a limited 
range, for LO between, say, 5 to 10 fm -1 . Variation of 
the cutoff in this limited range would not exhibit any of 
the drawbacks pointed out here, because it leads to an 
effective suppression of the higher-order terms. There- 
fore, reasonable fits to the data can be achieved. This 
explains the success of existing fits 0, 0, 0, 0, |2(j 
over a limited cutoff range. 

The variation of the cutoff within a given range will 
generate a band of values for observables. The error in a 
fit based on Weinberg's counting is likely dominated by 
the lowest partial wave without the required counterterm. 
As one goes to higher orders in Weinberg's counting, one 
acquires more counterterms, pushing the error to higher 
waves. The I suppression then ensures that the bands for 
the observables shrink, as observed 0, 0, 0, 0, |2(| • 
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FIG. 16: Comparison of attractive triplet phase shifts (as function of the laboratory energy) for A = 20 fm 1 (solid line) to 
the Nijmegen PWA (dashed line). 
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FIG. 17: Cutoff dependence of the triton binding energy. 

This does not, however, imply that Weinberg's power 
counting is correct. 

At any given order, our modified power counting has 
more short-range parameters than the same order in 
Weinberg's power counting. In this context, it is inter- 
esting to note that existing results in Weinberg's power 
counting already suggest that short-range parameters are 
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FIG. 18: Numerical results (crosses) and fit (solid line) for 
the deviation of the triton binding energy from its converged 
value, as function of the cutoff. See text for details. 



relatively more important than the long-range physics 
that is supposed to be of the same order. For exam- 
ple, Refs. ]r^.ll8| considered N 2 LO in Weinberg's power 
counting plus D-wave counterterms and found good fits. 
These calculations have a couple of counterterms more 
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than what we advocate here, but they come pretty close 
to our N 2 LO. 

A similar observation can be made about the Nijme- 
gen PWA when we look at it from the point of view of 
chiral EFT, with the PWA short-range parameters play- 
ing the role of counterterms. The long-range strong in- 
teraction consists of OPE and (leading and subleading) 
TPE [H H3, and thus corresponds to N 2 LO in Wein- 
berg's power counting. In the Nijmegen PWA, short- 
range parameters are added for the various partial waves 
until the fit to the data, up to 350 MeV laboratory en- 
ergy, is optimal. (The point is nicely illustrated by Fig. 1 
of Ref. [24(, where the quality of the 3 Po phase shift is 
shown for an increasing number of short-range parame- 
ters; see also Figs. 2 and 3.) The number of "counter- 
terms" needed per partial wave is larger, however, than 
in Weinberg's power counting in N 2 LO. In fact, if we as- 
sume that some parameters are needed for "fine tuning" 
to the data (c/. again Fig. 1 of Ref. [HI), or that they 
would not be needed if one were fit to the data only up 
to a lower energy (say, 250 MeV), then it appears that 
the number per partial wave is closer to what Weinberg's 
power counting would prescribe in N 3 LO. It would be 
very interesting to make this analogy between the Nij- 
megen PWA and chiral EFT more precise by fitting the 
data with a number of counterterms mandated by the 
different power-counting schemes. 

Another related point, which deserves further atten- 
tion, is that a good description of the AW data was ob- 
tained by the N 3 LO interactions in Weinberg's counting 
l2fj| . These calculations automatically include the 
counterterms in all partial waves that we consider to be 
non-perturbative. It is encouraging to see that a good de- 
scription of the data is obtained. However, lower cutoffs 
were employed than our LO study would suggest. This 
surprising fact can be understood if the range of cutoffs 
for which converged results are obtained increases toward 
lower A in higher orders of the expansion. This would be 
consistent with the observation that the triton binding 
energy is very well described for small cutoffs in Wein- 
berg's NLO. This needs to be studied more carefully in 
the future. 

We have taken here the minimalist point of view that 
only counterterms that are infrared enhanced with re- 
spect to naive dimensional analysis should be promoted. 
Since we have found no significant cutoff dependence in 
the 3N system in LO, 3 A forces are not infrared en- 
hanced at this level, and could be considered subleading. 
The same is true of the effective range in the 1 So NN 
channel. For both the 1 So NN phase shift and the tri- 
ton binding energy, our results are internally consistent. 
However, they are also somewhat unsatisfactory when 
compared to experiment. Agreement should be improved 
in subleading orders. An alternative, less-conservative 
approach would be to invoke a promotion of one or both 
of these interactions to LO on the basis of fine tuning. 
The relatively large value of the effective range supports 
this viewpoint. However, this is largely an issue of con- 



venience that we leave to later investigation. 

The LO results in our power counting are the ones 
given here. As we argued, there already exists some evi- 
dence that this power counting is consistent with pre vious 
results. A calculation beyond LO is in progress j65j. 



VI. CONCLUSIONS 

In conclusion, we have reanalyzed the predictions of 
chiral perturbation theory in the AW system using Wein- 
berg's original power counting in LO. We have identified 
that the singularity of the tensor interaction is respon- 
sible for a significant cutoff dependence in partial waves 
where it acts attractively. Furthermore, we have shown 
that the addition of one counterterm in each of these 
partial waves removes this cutoff dependence. It also im- 
proves the description of the data, which we see as a 
confirmation of our approach. 

For the 3A binding energy, we found cutoff- 
independent results in the limit of large cutoffs. We con- 
clude that the finite range of the interaction prevents 
the Thomas collapse of the 3A bound state. Meaning- 
ful predictions for the 3A binding energy, which are not 
possible using pionless effective field theory, are possi- 
ble here. We emphasize that this does not invalidate or 
compromise the pionless EFT approach, but reflects the 
stronger physical constraints built into ChPT. 

Our approach and its relation to the previous work can 
be best exemplified in the 3 D2 partial wave. We found 
that the addition of a counterterm ensures cutoff inde- 
pendence and increases accuracy. But we also found that 
in some ranges of the cutoff the low-energy description is 
equally good without this additional counterterm. This 
reconciles the traditional approach with our new results. 
It seems nevertheless advisable to promote counterterms 
in some partial waves. This ensures that all partial waves 
are cutoff independent for the same A and improves the 
description of the data over a wide range of energies. 

We have discussed how these results can be understood 
from a power counting that includes angular-momentum 
suppression. This improved power counting suggests an 
ordering of interactions that is similar to Weinberg's 
power counting, except for the infrared enhancement of a 
few of the counterterms that contribute to lower partial 
waves. 

Our study clearly has to be followed by further inves- 
tigation in at least two directions. First, one would like 
to understand in more detail the interplay of scales in 
infrared-enhanced counterterms. A more detailed analy- 
sis of the renormalization-group running and limit-cycle- 
like behavior of these interactions could shed light on this 
issue. Second, one would like to carry out a similar inves- 
tigation for the new NLO and N 2 LO interactions. One 
would like to confirm that a good description of NN data 
can be obtained already at N 2 LO. Few-nucleon systems 
should also be reexamined, as the triton binding energy, 
for example, does not come out well in LO. We consider 
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these to be important remaining issues that need to be 
studied for a consistent understanding of the application 
of ChPT to nuclear systems. 
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APPENDIX A: PARTIAL- WAVE 
DECOMPOSITION OF OPE 

The operator form of the OPE potential is 
, r ( -* -*f\ 1 ( 9a\ 2 q-u x q-u 2 , . 

where q '= p ' — p' is the momentum transfer. 

The isospin operator separates easily. It is given by 



(t\\r 1 ■ r 2 \\t') = (2 t(t + I) - 3) 5 W 



(A2) 



The spin-orbital part can be decomposed into spin and 
orbital tensor operators using obvious basis states, start- 
ing with 

(p(ls)jm\ ^if^'f 2 \p'(l's')j'm') 
q -+- m^ 



I 00 



3 Ws)jm\ {gl Q l , {a l Q} \p'{l's')j'm') 



q 1 + m£ 



(A3) 



where we used the representation a-b = — y/Z {a b} 00 of 
the scalar product of two vectors a and b. With (q)\ = 

g Jj Yi\(q), we can recouple 



(p(ls)jm\ q :°' q '° 2 \p'(l's')j'm') 

i ^ _L_ <m ^ 




> (1 1 /,00) (j'O j,m' 0m) 



ls(pl\\ 



4irq 2 



q 2 + ml 



Yf(q)\\p' l') 



Ml {<T1 °2} f \\S>), 

thus separating spin and orbital parts. 
The orbital part is 

(plm\ Y fM I P' I' m' ) 



(A4) 



q l + m% 

I 47T/! 



Ai+A 2 =/ 



Ai! A 2 ! 



7 { 
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Ai A 2 / 
I V f 



> feVAiAz (k Ai 1,00) (k A 2 I', 00) 



(/' fl,m'^m) (-) 



I' .if 



(A5) 



where the angular dependence of the propagator was ex- 
panded in Legendre polynomials using 



9k 



f ( PP ') = ^(-) k 4n ^ dxPk{x) 



q 2 1 

Z7 T-(A6) 



q* + qJ 

This confirms that the orbital part is a tensor operator 
of rank /. 

With the spin matrix element 



} . (A7) 



f— \ 3 S ' f 

(s\\ {ai a 2 } } \\s') = 6yjs'f I 1/2 1/2 1 

[ 1/2 1/2 1 

one obtains the complete matrix elements of the OPE: 
(p(ls)jm\ q '° iq '° 2 \p'(l's')j'm') 



= E 



q 2 + ml 

1 1 / 

v47r I 



3V/ 
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I V f 
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APPENDIX B: REMOVAL OF SPURIOUS AW 
BOUND STATES 

To remove a spurious bound state in the AW system 
we change the interaction V to 



V = V + | X ) A ( X | , 



(Bl) 
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where \x) is the wave function of the spurious bound 
state, and A is an energy parameter, which determines a 
shift of the binding energy for the spurious state. The 
limit A — > oo removes the spurious state, and V is phase- 
shift equivalent to V. 

If t and i solve the LS equations for V and V, respec- 
tively, they are related by 

i=t+\ V ) N ( V \ , (B2) 

with 

\V) = \X) +tG \ X ) (B3) 

and 



Here Go is the free, two-particle Schrodinger propagator. 
This formulation allows us to perform the limit A — > oo 
analytically, and we end up with the i-matrix 



which then enters our calculations for the triton binding 
energy. The accuracy of this prescription can be checked 
numerically by comparison with the expectation value of 
the Hamiltonian H using the original potential V. 

This procedure can easily be generalized to two and 
more spurious bound states. 
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